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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10-3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

C. 01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

‘To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  =  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  =  (5/9)  (F  -32)  +  273.15. 
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ESTIMATION  OF  WAVE  TRANSMISSION  COEFFICIENTS 
FOR  PERMEABLE  BREAKWATERS 

fry 

Uillicun  N.  Seelig 

I.  INTRODUCTION 

The  primary  purpose  of  a  breakwater  is  to  reduce  wave  energy  in  an 
area  to  be  sheltered.  One  of  the  important  characteristics  of  a  break¬ 
water  is  the  magnitude  of  the  wave  transmission  coefficient,  defined  as 
the  ratio  of  the  transmitted  wave  height  to  the  incident  wave  height. 

Two  basic  types  of  wave  transmission  are:  (a)  by  overtopping  that 
occurs  when  wave  runup  exceeds  the  crest  elevation  of  the  breakwater, 
overtops  the  breakwater,  and  generates  waves  in  the  lee;  and  (b)  through 
a  permeable  structure  that  occurs  because  some  of  the  wave  energy  is  not 
dissipated  by  the  breakwater  and  is  transmitted  through  the  breakwater. 
The  total  wave  transmission  coefficient,  Ky,  is  given  by: 

Ky  =v/,(K  To)2  +  (Krt)2  =  ht/Ni  (1) 


where 

=  transmission  by  overtopping 

coefficient 

Kn 

=  coefficient  of  transmission 

through  the  structure 

Hr 

J. 

=  incident  wave  height 

Hy 

=  transmitted  wave  height 

These 

and  other  symbols  are  defined 

in  Figure  1. 

Since  the  prediction  method  is  complex,  particularly  for  transmission 
through  the  structure,  a  computer  program  is  presented  in  an  Appendix  to 
this  report.  The  program  incorporates  the  analytical  model  to  determine 
Ky£  by  Madsen  and  White  (1976) 5  and  an  empirical  equation  to  determine 
Ky0  developed  by  Seelig  (in  preparation,  1979)2. 


iMADSEN,  O.S.,  and  WHITE,  S.M.,  "Reflection  and  Transmission  Character¬ 
istics  of  Porous  Rubble-Mound  Breakwaters,"  MR  76-5,  U.S.  Army,  Corps  of 
Engineers,  Coastal  Engineering  Research  Center,  Fort  Belvoir,  Va.,  Mar. 
1976. 

2SEELIG,  W.N.,  "Two-Dimensional  Tests  of  Wave  Transmission  and  Reflec¬ 
tion  of  Laboratory  Breakwaters,"  U.S.  Army,  Corps  of  Engineers.  Coastal 
Engineering  Research  Center,  Fort  Belvoir,  Va.  (in  preparation,  1979). 
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Ky  =,/(Ky0)  2  +  (Kyt)2  -  Ily/llj 


Figure  1.  Definition  of  terms  for  wave  transmission 
for  permeable  breakwaters. 

II.  WAVE  TRANSMISSION  BY  OVERTOPPING 


Wave  transmission  by  overtopping  occurs  when  wave  energy  is  trans¬ 
mitted  by  flow  over  the  top  of  a  structure.  The  transmission  by  over¬ 
topping  coefficient  can  be  estimated  using  (Seelig,  in  preparation, 
1979) '  : 


KTq  =  C  (1-F/R) 

=  0  for  F/R  greater  than  1.0 


(2) 


where 

R  =  wave  runup 

F  =  breakwater  freeboard,  defined  as  the  structure 

height,  h,  minus  the  water  depth,  d„ 

1  u 

C  =  an  empirical  coefficient 


(Ky0)max  =  1.0 

Laboratory  tests  show  that  the  value  of  C  is  related  to  the  crest  width 
of  the  structure,  B: 

C  =  0.51  -  0.11  B/h  .  (3) 


Thus,  a  slight  decrease  in  the  transmission  by  overtopping  occurs  as  the 
structure  crest  width  increases. 


3SEELIG,  IV. N.,  op.  cit.,  p.  7. 
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Wave  runup  is  estimated  using  the  formula  (Ahrens  and  McCartney, 
1975) 4 : 


R  _  a£ 
%  1  +  b£ 


(4) 


where  a  =  0.692  and  b  =  0.504  are  recommended  for  rubble -mound  break¬ 
waters  and  a  =  0.988  and  b  =  0.703  are  recommended  for  a  breakwater 
armored  with  two  layers  of  dolos.  £  is  the  surf  parameter  given  by 


where  9  is  the  angle  of  the  seaward  face  of  the  breakwater,  and  L0  is 
the  deepwater  wavelength  obtained  from  linear  wave  theory.  Calculations 
of  wave  transmission  by  overtopping  are  performed  automatically  in  the 
program  MADSEN  (see  App.). 

III.  WAVE  TRANSMISSION  THROUGH  PERMEABLE  BREAKWATERS 

The  coefficient  of  wave  transmission  through  permeable  breakwaters, 

K Tt>  is  estimated  using  tne  analytical  model  of  Madsen  and  White 
(1976)5.  in  this  model  the  transmission  coefficient  is  related  to  a 
complex  function  of  the  size  and  porosity  of  the  materials  used  in 
building  the  breakwater  (Table  1),  the  breakwater  geometry,  the  seaward 

Table  1.  Porosity  tf  various  armor  units  (after  U.S. 


Army,  Corps  of  Engineers,  Coastal  Engineer 
ing  Research  Center,  1977)6. _ 


Armor  unit 

No.  of 
layers 

Placement 

Porosit.v 

Quarrys tone (smooth) 

2 

random 

0.38 

Quarrystone (rough) 

2 

random 

0.37 

Qua  rrystone (rough) 

>3 

random 

0.40 

Cube (modified) 

2 

random 

0.47 

Tetrapod 

2 

random 

0.50 

Quadripod 

2 

random 

0.49 

Hexapod 

2 

random 

0.47 

Tribar 

2 

random 

0.54 

Dolos 

2 

random 

0.63 

Tribar 

1 

uniform 

0.47 

Quarrystone 

graded 

random 

0.37 

‘‘AHRENS,  J.,  and  McCARTNEY,  B.L.,  "Wave  Period  Effect  on  the  Stability  of 
Riprap,"  Proceedings  of  Civil  Engineering  in  the  Oceans/III ,  June  1975,  pp. 
1019-1034  (also  Reprint  76-2,  U.S.  Army,  Corps  of  Engineers,  Coastal  Engi¬ 
neering  Research  Center,  Fort  Belvoir,  Va. ,  June  1976,  NTIS  A029  739. 

5 MADSEN,  O.S.,  and  WHITE,  S.M.,  op.  cit.,  p.  7. 

flU.S.  ARMY,  CORPS  OF  ENGINEERS,  COASTAL  ENGINEERING  RESEARCH  CENTER, 

Chore  Protection  Manual ,  3d  ed.,  Vols.  I,  II,  and  III,  Stock  No.  008-022- 
00077-1,  U.S.  Government  Printing  Office,  Washington,  D.C.,  1977. 
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slope  of  the  stricture,  water  depth,  wave  height  and  period,  and  the 
kinematic  viscosity  of  water  (Table  2).  To  use  this  method,  waves 
should  have 

d  1,25 

L  VHj  cot2  0 

where  L  is  local  wavelength. 


Table  2.  Kinematic  viscosity  of  water. 


Water  temperature 
(°C) 

Kinematic  viscosity  of  water 
(m2/ s) 

0° 

0.0000018 

10° 

0.0000013 

20° 

0.0000010 

30° 

0.0000008 

The  Madsen  and  White  model  was  tested  against  laboratory  data  for 
permeable  breakwaters  (Seelig,  in  preparation,  1979)7  and  was  shown  to 
give  useful  estimates  for  both  monochromatic  and  irregular  waves.  For 
irregular  wave  conditions,  the  wave  input  to  the  program  should  be  the 
mean  wave  height  and  period  of  peak  energy  density.  A  few  tests  with 
breaking  waves  suggest  that  the  prediction  method  can  also  be  used  with 
breaking  waves.  The  Madsen  and  White  model  appears  to  effectively 
account  for  breaking  wave  energy  losses,  although  it  does  not  explicitly 
include  breaking.  Tests  of  breakwaters  armored  with  dolos  units  suggest 
that  the  program  can  also  be  used  for  artificial  armor  units.  Compari¬ 
son  with  laboratory  data  shows  that  the  model  gives  the  best  predictions 
for  shallow-water  waves.  Predictions  of  transmission  coefficients  tend 
to  be  conservative  for  transitional  or  deepwater  waves.  Refer  to  Seelig 
(in  preparation,  1979)7  or  Madsen  and  White  (1976)8  for  more  information. 
Figure  2  shows  a  comparison  between  wave  transmission  coefficients  ob¬ 
served  in  a  laboratory  model  and  predicted  using  the  methods  described 
in  this  CETA. 


IV.  EXAMPLE 

Use  of  the  computer  program  (MADSEN)  in  the  Appendix  can  best  be 
illustrated  by  an  example.  The  format  of  required  input  information  is 
given  in  Table  3.  Any  number  of  breakwater  geometries,  water  depths  or 
wave  conditions  can  be  analyzed  in  a  single  run.  The  first  53  cards  are 
a  standard  deck  of  look-up  (input)  tables  (see  Table  A- 1) ;  card  type  1 
provides  the  number  of  breakwater  configurations  or  water  depths  to 
analyze.  Card  types  2  to  6  give  required  input  information  for  each 
breakwater  of  interest;  however,  a  separate  set  of  these  card  types  is 
required  when  the  breakwater  geometry  or  water  depth  is  changed. 


7SEELIG,  W.N.,  op.  cit.,  p.  7. 

8MADSEN,  O.S.,  and  WHITE,  S.M.,  op.  cit.,  p.  7. 
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Table  3.  Format  of  input  information . 


Card  type 

Format 

Description 

Standard 

53  standard  input  cards 

1 

12 

Number  of  breakwater  configura¬ 
tions  or  water  depths  to  test 

2 

20A4 

Title  card 

3 

312, 4X, 7F10. 5 

Number  of  wave  conditions  to  test 

Number  of  materials 

Number  of  horizontal  layers 

Structure  height  (m) 

Water  depth  (m) 

Kinematic  viscosity  (m2/s) 

Width  of  top  of  breakwater  (m) 

Front  slope  of  breakwater  =  tan 

C©) 

Wave  runup  parameter  a  =  0.692 

Wave  runup  parameter  b  =  0.504 

4 

lOx, 2F10. S 

(1  card  per  material) 

Material  diameter  (m)  (armor  1st) 

Material  porosity 

5 

10x,7F10.5 
(1  card  per  hori¬ 
zontal  layer) 

Layer  thickness  (m) 

Mean  length  of  each  material  type 
in  the  layer  (put  in  consecutive 
order;  e.g.,  material  1  (armor) 

1st,  etc.) 

6 

2F10.5 

(wave  condition  card; 
one  card  per  wave 
condition) 

Wave  period  (s) 

Wave  heights  (m) 

Repeat  card  types  2  to  6  for  each  water  depth  or  breakwater  configura¬ 
tion  to  be  tested. 


Card  type  3  gives  the  number  of  wave  conditions  to  analyze  and  sum¬ 
marizes  general  input  information  (Table  3).  For  the  example  breakwater 
(Fig.  3),  18  wave  conditions  with  periods  of  5,  10,  and  20  seconds  and 
with  heights  that  range  from  0.1  to  2.0  meters,  are  analyzed. 

Card  type  4  gives  material  characteristics,  one  card  per  material 
and  the  first  card  should  describe  the  armor  material.  The  example 
gives  three  materials  (armor,  underlayer,  and  core);  diameter  and  poros¬ 
ity  of  the  materials  are  shown  in  Figure  3. 

Card  type  5  is  used  to  input  the  mean  horizontal  length  of  various 
materials  in  various  horizontal  layers  of  the  breakwater.  A  new  hori¬ 
zontal  layer  occurs  when  there  is  a  change  vertically  in  material  type 
or  slope  and  the  layer  next  to  the  seabed  should  be  designated  as  "layer 
number  1."  In  the  case  of  the  example  breakwater,  three  horizontal  layers 


12 


are  shown  in  Figure  4.  Sample  horizontal  length  calculations  are  also 
included.  Note  that  when  determining  horizontal  lengths  of  the  armor 
material,  the  outer  layer  of  the  armor  on  the  seaward  side  of  the  break¬ 
water  should  be  "removed"  first  because  dissipation  of  the  seaward  face 
is  determined  in  a  separate  part  of  the  computer  program. 


Table  4  gives  the  computer  program  input  information  required  for  the 
example;  Table  5  is  the  resulting  program  output.  The  output  shows  that 


predicted  transmitted  wave  height  for 
of  incident  wave  height  and  period. 

Table  4.  Sampl 

1 


ex  »«pi  e 

PROhUH 

18  J  J 

6,0 

«.& 

000O009S 

WAT  1 

0.729 

0.  IT 

WAT  2 

0.339 

C.J7 

HAT  J 

0.0<>2 

0.37 

LAV  1 

3.55 

6.51 

3.80 

LAV  2 

0.7  A 

6.53 

2. 56 

LAV  1 

o.uT 

5.25 

0.0 

5.0 

fl.l 

5.0 

0.5 

5.0 

1 .0 

5.3 

1.5 

5.0 

1.75 

5.6 

2.0 

10.0 

0.1 

10.6 

0.5 

10.0 

1  .0 

10.0 

1  .5 

1  o  ,  0 

t  . 75 

lo.o 

2.0 

20.0 

0.1 

20.0 

0.5 

20.0 

1.0 

20.0 

1  .5 

20.0 

1  •  75 

20.0 

2.0 

this  example  is  a  complex  function 


input . 


2.52  0,667  0.692  0.50 a 


6,00 

0.0 

0.0 


V.  SUMMARY 

A  computer  program  is  presented  for  estimation  of  wave  transmission 
coefficients  for  permeable  breakwaters.  Extensive  testing  of  the  program 
with  laboratory  data  has  shown  that  the  program  can  be  used  to  estimate 
transmission  coefficients  for  monochromatic  or  irregular  waves  and  for 
rubble-mound  or  other  types  of  permeable  breakwaters.  A  limited  amount 
of  testing  suggests  that  it  can  also  be  used  for  breaking  and  nonbreaking 
waves . 


A  copy  of  the  card  deck  and  more  extensive  program  documentation  for 
the  computer  program  MADSEN  (CERC  Program  Number  752X6R1CPO)  are  avail¬ 
able  from  the  ADP  Coordinator  at  CERC.  The  cost  of  running  the  program 
on  a  CI)C  6600  computer  is  only  a  few  cents  for  each  wave  condition  of 
interest . 
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Table  5.  Sample  output. 


Example  PROBLEM 

COMPUTATIONS  Of  wAVf  TRANSMISSION  THROUGH  A  POROUS  SRCAKNATfR 


NUH  OP  WAVE  CONDITIONS  18 

NUM  OP  MATERIALS*  3 

NUM  op  hOHUONTIAu  LAYERS*  s 

STRUCTURE  HEIGHT  <m>«  6,60* 

WATER  depth  fH)s  U, *00 

WINEMATIC  VISCOSITY  (M2/slC>»  ,000000630 
Sw  TOP  WIDTH  (m)j  2.520 

TANS  op  PROMT  slope*  .6670 


runup  COEFFICIENTS  A*  .602  6*  ,5flO 
material  characteristics  (m*ke  armor 


material  number  n 


MATERIAL*  1  OIAMfTER  (mj«  ,729  POROSITY*  .JTO 
MATERIAL*  2  OTAMETtH  I *  .33*  P08D8ITV*  *  3  7  0 
MATERIAL*  i  01 AMfc TfcR  (M)«  ,002  POROSITY*  .370 


HORIZONTAL  LAYER  CHARACTERISTICS 

(make  LAYtR  NEXT  TO  SEABPU  LAYER  NUMKER  1) 


MATERIAL* 


HORIZONTAL  layer*  1  THtChnPSS  (Ml*  1.550  lengths  (M)» 
MORIZONTIAL  LAYE«*  2  ThtChnESS  (M),  ,?*0  LENGTHS  (Mj, 
HORIZONTAL  LAYER*  5  THTCHNESS  (Ml.  .070  LENGTHS  (H), 


1 

a.* 
a. 5 
5.3 


2 

1.* 

2.5 

0.0 


1 

6. a 
0.0 
0.0 


H(H1 

T(SEC) 

H/(G*T*f) 

.too 

5,00 

.OOOOflO 

.500 

5.00 

,002001 

I.ooo 

6.00 

,0000*2 

1,500 

6.00 

,006122 

1.750 

5.00 

,0071a J 

2.000 

6.00 

.00*16 j 

.100 

10,00 

.000102 

.500 

10.00 

,000510 

1.000 

10.00 

,001020 

1.500 

10,00 

.001531 

1.750 

10.00 

,001 7*6 

2,000 

10.00 

.002001 

.100 

20.00 

.000026 

.500 

2o,  00 

,000126 

1.000 

20,00 

,0002*5 

1,500 

20,00 

,0003*3 

1.750 

20.00 

•  00  0  006 

2.000 

20.00 

.ooosio 

N/L 

D/(G*T*T/) 

kit 

ktd 

KT 

,00335 

.01 66 

.3*1 

0.0*0 

.391 

,0l67o 

•  0  1 66 

.211 

0.000 

.211 

.03309 

.0196 

.109 

0.000 

.109 

.05023 

,01  *6 

.126 

.036 

.130 

.05660 

,0166 

.121 

.0*6 

.106 

,06697 

,01*6 

.113 

.125 

.166 

.00151 

,00«9 

.367 

0.000 

.397 

.00753 

,0009 

,199 

0.000 

.169 

.01507 

,0009 

.135 

0.000 

.135 

.02260 

,00«6 

.099 

.115 

.02637 

,0009 

,oaa 

.169 

.1*2 

.03013 

,0009 

.0*0 

.163 

.209 

.00073 

.0012 

.379 

0.000 

.376 

.00367 

.0012 

.1*0 

0.000 

.1*0 

,00735 

.0012 

.125 

.010 

.l** 

.01102 

,0012 

,0*6 

.150 

.1®* 

.012*6 

.0012 

,0*6 

.196 

.216 

.01070 

.0012 

.060 

.227 

.201 

Kfi  MT ( Ml 
,26  .030 
.28  .105 
,26  . 1 ao 
,27  .201 
,27  .200 
.26  .337 
.51  ,oao 
,60  ,000 
,62  .135 
.60  .228 
,60  .316 
,60  ,01* 
.53  ,036 
.66  .002 
.70  ,125 
.71  .273 
,72  .376 
,72  ,061 


KTt  •  wave  transmission  Through  the  structure 

KTo  •  WAVE  TRANSMISSION  «y  nVfRT(jRPlNG  COEFFICIENT 
KT  •  TOTAL  wave  TRANSMISSION  COEFFICIENT 
KR  .  WAVE  REFLECTION  COtPPlCIENT 
NT  •  TRANSMITTED  WAVE  Mt I GhT 


■lit  MU  IS  —  •  •  il  i 

|UOM  9<tf I  *v  ,;uC  — *’*'*" 
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APPENDIX 


LISTING  OF  THE  COMPUTER  PROGRAM  MADSEN 


Eo« 

PHOG«*m  MAr>Sf  M(  ImPiiT.ijiiTPuI  ,  TaPMsTNPUT,  T  APk6»0UTPUT  ,  TAPk  51 

CUMMON/MAUSI  /NM,ML.i'(  1 1  )  ,Mf  j  t )  ,Ll  ( 1  J  <  1 1  )  .  TH(  1 1 ) 

coMKUN/sm  /nkl.fs 

RkAL  M«l 

f)JM|  MSICN  IHMF  (  J  )  ,  T  I  TLF  (PO  )  .  NiiM(  1  0) 

RML  IiM'*HI,KR,kiiU  •nR*LL  .KTO.KTT 
OAT  A  M.lM/  S  ,  ?  ,  1  •  U  ,  b  >  b  .  7,0,9,  1  0/ 

P)  =  3.  M1S9 
CALL  RtAiil 
RMDt'xSOO)  wCO'-p 
590  FORMAT  f  i  T  ^  .  ■m,  7r  t  O.S) 

On  2f>0  1 .1=  1  .  scO-ip 
C  RfAo  input  information 

R|.  AO(5*  I  71)  f  T  I  n  1 1  J.JM)  ,?(>) 

1  7  1  FORMAT  C  ^0  Ail  ) 

HRITtCF.fi/?)  fTlTLF(JjM),JJM3t,20) 

J7j  FORma)( JH| , |/>x.?0A<i) 

RMt)(SfS«»1  NTf  NM,NL,HS,MUfMi,  T{)PW,  T  ANR.RA.RB 

f  S  H  5)  •  H  ( • 

TFfRA.ll.O,)  WAsn.AV? 

IMRM.It.O.)  RBs.SOu 

wRTTf  (hi  *11 1  )  NT  ,  \M,  Nl  «  HS,HO,  On,  TOP-.  T  AOh«RA, 

97l  FORMAT  f  /  ,!Oa.  'rU-PuTAllORs  Of  RAvl  TrAK>>"ISSION  THROUGH  A  POrOiiS 

♦  HR)  «t>.  »  a  1 1  a  i  , /// ,  s«  i  »mii«  pF  ..  a  vF  COwOM  IONS' •  t^X,  1  3./.SX, 

♦  IM/M  (IF  MATFR1AI  SS'. 17A.M. /.Sx, 

♦  inijm  Of  «()RT^(/ Vt;  At  i  Avfa-Ss  *  ,  A,x.  ISf/.SX.  'STRuCTURfc  MfTGHT  Cm) 
»S*.6X,»  10.3./.S*.  'AATtR  PfPTM  O  Is  >  ,  I  1  X.F  1  0,  $, /.*>X. 

AtKIK.Fl.ATlC  vtsCosITt  (M?/bf  C  )=  I  ,f  1 1  ,R,  /,SX,i«K  T  (IP  w  ]  0  T  u  (M)S1, 

♦  10X.FI  o.  3. /,sx.  '  TA'<0  Or  FRONT  S|.(1P£S '.  9X.F  «,«,/,  Sx  RUNUP  ct/fFr  Icl 

♦  FhTS  »s'.rb,J«'  «3'«kb.3) 

00  99  I  =  1»U 
DO  VI  .1=  1  .  1  1 
90  U(I»J)  =  0, 

99  CONTlNlJt 

WRIT!  (K..?«3) 

*«3  FORMAT  ('j*.  imaTFRIAL  C  *  A*  AC  T  k  R  I  ST  T  C  S  (MaKF  ARMOR  MaURJAI  MlMHkR  l) 

•'«./) 

00  b  1 = 1 . N» 

Rfc  ADC  S • 7 1  n c  I ) .  n  f  I  > 

7  FC.wt'ATCIfx,  7Fift,«.) 

WKMt.C6.177)  I.ncM.NCT) 

177  FowmAtCSX,  tMATF  ol  At.e  •  ,  Ti,  *  OTA-tTfR  Cm)*i,F  6,3.'  POROSIT  Ys  '  .  t  b  •  3 ) 
b  CONI J "OF 

WKJTtC6.^»«)  C  Ntj*'C  JM)  .  f>  s  1 

2«a  FORMATC//.SX.  *M<  RI/uMtai  l  AVf«  r  H  A  w  A  C  T  k  R  J  S  T  I  C  S  1  I  /  .  S  X  , 

♦  •  C  MA*f  1  a  r  F  R  M  XT  To  SrAMfO  LavIR  NOMRFW  1)'f/. 

♦  5?X. »maTfbial»  ’»7C  |1 .Sx) . / . 6 3 x . 6 ( I?. u X ) ,/) 

00  $\  J* 1 .nl 

Hk  AO  C  S •  /  )  T  t.  C  J  )  ,  M.  L  C  T  .  .! )  «  J  e  J  ,  k,m  ) 

♦R  I  Tk  Cfr. t  7A)  J.TwCJ).(lLCl.J).’8l,NR) 

17«  FOWATCS*.  imOrMonTUl  LAYF««t,l3,t  TnicmnfsS  (Ml*tt  F6.3.I  LkAGTM 


gu  rm  u  nsi  ?wnx  yt  nuss-i  *  * 
ntoM  oopy  ?»■  .or. 


*S  fM)sit7F*.tt/'ht>Xt7tfc,i) 

55  CONTINUE 
NMSNmA 1 

n<M*)ar>{  t ) 

N(  N*) eO.it  I 
NLsnLa 1 

TH(»4L)sjnoMiftoP# 

» =  5 .  •  f  t  ) 

9«2  FURM*TF//.  bX,'NF*)  T(Sft)  H/(r,*T*T)  H/L  r>/(C*T»r/)  KtT 

«  K  T  «>  M  K»  Hi  f N)  I  ) 

On  l 9g  IKsi.nT 
RLAnFs*«>  T,m 

e  ruRMit to.s) 

A  8  H  ♦  0  ,  S 
OHsD( 1 )*0,S 

IF  (*. I  T, 0,o » ft 01  )  GO  TO  mo 

IMTANh.U  .0.)  f,n  TO  37 

CALL  PFFL(A,MS.nFl>»M0,TANR,T,»I!fRU.L) 

AIs«II»A 

2?  DHTs?,*Mj*A 
IF  L  Al<:D 

C  *sSUML  l)*H  s  OH  t  AMI)  ITtRATt  ON  i  me  FfJliJIlVANT  Bw 

ICOilNTsy 
DHF  sl)t<  r 

10  1 1  OI)N  rs  I  COiinT  ♦  1 

CALL  {  OD"(OhI  ,  DhT  •  l  t  iHdi  M$f  T  AnB*  NR  ■  l)R.  TOP*  ) 

CAl  L  I  n  T  t  R  (  JR|  T  ,  lL«hii,  A  1  ,N(i,n«.  T  1  ,RI  f  L«  IF  L*0) 

I F  ( I F  l  a  (. ,  M; ,  1  j  npsnwAO.NS 
I F  ( l  F'  |  ad.eo.i)  <;n  TO 
oHfc- ( )  ,*wn»Ri  i*» 

JFdLOUNT.il.u)  GO  10  10 
KF<BI<UF<1 1 
KITst J*RI1 

37  IF  (  T  AnB.LF  ,o.  )  CALL  I N  U  R  (  N  (  1  )  ,  T  ,  T  OF'*.  ,  HO,  A  «  *0,  f)  (  |)  ,  K  T  T  ,  K»,  l  ,  I F  l  AG) 
IF  (  IFLAO.F.O,  1  )  DRSI)K*0.S 
Id  UL*O.FR.  1  )  GO  Til  3  7 
SuRF  =  T*NH/S«RT(H/( 1 , So* 1  *  T  )  ) 

RHSRA»51)HF/(1  ,  ♦R«*SURF  1 

Rs  H* Rw 
FRsF/R 

CeO.Si  -0.1t*TOF»rt/HS 
KT0sC*( I  ,-F  R) 

IF  F  (  ToF»*/»S)  ,(,T  .O.HH.AnO.f  ,LT,0,1  * T Use* < l ,-F R ) - F 1  * 1 1  * F « 

IF FkTo.GT.i  ,)  KTp=1, 

IF  FF«,01  ,  J  ,0  1  KTfisn, 

HC.T2Sa*2./(<f,ro*t*T) 

Ml S^,AA/L 

OGT^shO/FR.moaTatJ 
FLAGs  in 

KT«SRRT(kITa*?*kTU**2) 

IMKT.GT.t.O)  KTsl.o 
HI BH* A  T 

MRlTEFO.Rfll  )  M,  t,nGT2,mL  ,or,T?,t'TT,«in»KT,KRf  Mt 
98l  FoRnATF  SX,Fh,3*F  1  0,?«F  lO.ft.F  j  o,s,r  10,<jf  3F  6, 3,F6.?,F  7,31 
ton  continue 

199  CONTINUE 

MRITF  FO*?Ol 1 


an  <..■  -,«t  *«u>m  18 

now »n»  ••  -  '<riC 


Z0\  FORMATC//.?*.  IKtT  •  »AVL  TRANSMISSION  THROUGH  IMF  S  T  BUT.  T  URf  •  »  /  t 
*?X  t  •  XTO  *  w*VF  T  R  A  N  S  M  1  S  S  T  0  '■/  H  Y  0V{.  R  1  t'PP  I  h(.  C  Ut  F  F  !  C  1 f  N  T  1  «  /  • 

*  i  l  k  T  <■  TOTAL  «*vF  ThanSmTsSJOn  CULM  ic  U  Nf  i  , /,?», 

*  •  k w  -  »avf  HtPifct ton  cofff jr if nt  •  • 

*  /  » <?X  #  I  m  •  TWANSHlTUn  »AVt  hFI(.MM) 

200  CO^TIMifc 

STOP 

F*0 

SURRCHlTIU  RFFI.  (A.Hb.l'.HO,  TANg.T  ,BT  I.  hi  l,  l  ) 

COH^OH/HAUS/KST (OiII)i rut (O, t)1,RTfl?#in«Tx(R,JO),KX(R|10T 
niMfA-sH-A.  F5R(  1  n  »HOS(  1  J  )  .RSTl  1  ) 

Rt AL  I  •  l  SL • L S 

C  CF  s  H0»M.  CoRRFCTIon  F  AC T UR  TCI  ACCOUNT  FOR  fiODF  L  SI  OPT  [FFTCtS 
CPs  1  ,?«-0,S7P»TakjH 

IKTAf.H.t  1  fTsl.O? 

IFtTAhH.Gt.o.hC)  CTso.aR 
C  r  JMD  WAVF  If  'Jl.tM  l. 

MOLOsmO/ (J.SH*T*T) 

CALL  LFMiTfHOLP«H»  l  ) 

Lsho/hUl 
LSshQ/I  a  nh 

IF(HS,H,ho)  iSshS/!Am« 

LSLsLS/L 

U  (1  S(  ,LT  GO  TO  US 

lMtNaSrTRT(H,?M*(t3/fl.4)/(<).A»TAFiH(t>t?l«J*HP/(lS/0(A)))) 
mRITF.  (b»  lot  )  T  M I  h 

lot  FORHATf///,lXt  '..arm1n(;.Tm£  m  I  a,  T  mim  »AvF  PI  R  1 OD  TO  HI  ANALYTIC)  R  V  T 

*  MIS  “w'H.mah  T  S  1  i  f  ( •  i'  •  *  Sft  FOR  This  CONUMIU*") 

L  SL  s  0 , 7  9R 

1  Os  Is(LSL*|o,M  ,) 

C  I  N  T  F  RPOL  ATI  I  •♦PUT  Ta«LL  FOR  THIS  LSL  ViLUfc 
Il»LSL*20,*l. 

00  i  J=t  » 1 t 

F$S< J)=F  ST  ( I. J)»(F  S  T  (  1*1 .  J)-TRT(  T  ,.)))♦  f  LSL- (  1-1  )*0,  l)/0,  J 

pus (  j i shut ( i .sn*(RiiT c  uu.n-RuT ( i,.n )a(lsl-u-i  >*o.n  /o.i 
J  RSfJTsRT  ( 1  T ,  J)«  (R>  fl  !♦  |  *  J)-RT(  TI#  J)  )*(LSt  -  U  UiUo.OS)  /0,0«T 
C  Gill-  SS  PHI  A.\0  ITFRAIF 
PmJsS.O 
Ms  n 

6  JaPHl 

FACsCaI  or,(PMlfl  .T-ALOOf  Jfi.)  )/(AI  OGC  Ja?.)-alOc;cJa|  ,) ) 

FS=F  SS(  J*  !)♦  F  AC*  {FSSC.I+?)-FSSC  J«  1  )  } 

RU=RUS(  0*  1  )  ♦  F  AC*(R|'S(  .!♦  ?  J -«« IS <  J4  1  )  ) 

RlIsRSf  J»n»CRS(.T*2)-RS(.lM  )TAF  AC 

ARr.s0,?V»(O/Mii)«»l).«J«(R''A2,*A/(Hn*TAMA)J»#0,i#rs 
PH  I  A.S0,S»  A  T  AH  {*RG)*S  7 
HSMf  J 

OtLsARSCPHtN.PMJ ) 

IF  (M.f.T  ,?01  GO  TO  R 
PH I  a  Pm  Is 

IF  CP“1  ,L  T.O.Ot  )  PHlso.m 
IF  (PHJ  ,0T  ,0,09)  PMlsR.OV 
IF (Ofc l ,GT .a.ftS)  GO  TO 
9  RIIsRlMCF 
ML  TcjRn 
FNO 

SuhrOuTIni  rfaoi 

COHHOn/mAOS/FSTCO,  t  n»RUT(R,1  1  *IX(9.I0)  .«X(R,10) 
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ai(  uu  i»  bcsi  «MAwn  nw**** 

nod  CO*  T  w  DG8 


17? 

1 

2 

3 

4 

b 


1 


3 
b 

4 


>Oo 

2 


21 


4 


5 


J7 


f 0«M*t (  S<«  fTFtt.iM 

no  t 

»M0(*>*17/)  ffSl(N.M)«N*l«9) 
DO  2  M*1 *  1  1 

HMD  (  S  «  1  /  7  )  f  Wilt  f  N.M)  ,w»1 ,9) 

no  3  m«1 •  1 1 


Bf  AD (5077)  («T(M»m)iNs1*1?) 

no  4  MS  1  ,  1  0 

Rt  AO ( 6  O  7  7 )  {TX(MtM)*i«sl,9) 

no  s  m=  i ,  i 

«tAO(«>#i77)  (k*(n."i #*«sr.<n 

»t !o«0 

{•NO 

sumwooT 1st  t g ng i r 
MtAL  l  !•»  I  O^L*  *  l  (111 
IP=1 .0/01 0 
LOPS  1  («/iHn 
NS  1 


Pl*1.1“ts« 

AM6s^.ii*P  |  /lo 
1,0^1  *sLon»  T  ANh(  APG) 

NS  N  ♦  t 

OIFF  =  AMS(LONF*»l.n) 

IX(N-i>O0)  I, <1,4 

If  (MO -o.oons)  ? * ? * 5 
L0=(L0N1  -♦1.01/2,0 
GO  10  1 

m  =  i  .o/i  o'*e« 

*4IH|hiioA)  run, in. 

FuWf'4  T  (  41|M  SnuMOiillNt  I  t  m;Tm  0 1  0  Mil  COMVfHGti  D/LO 
1  OHO/i  s  , F 1 0 . b ) 

nt  =  i  .(i/u>Nf* 

T U w n 
F  NO 

sopwoo 1 1 Mt  f of»w( not . omt. 1 1 «nn, hs. r anh,n«* l>w. To*’* ) 
COMMON /"A  OS  1  /NM,  Ml.  .0(1  I  )  .  N(  J  1  )  .  L  C  1  1  I  1  1  )  .  1  H(  1  1  ) 

0  I  Mf  NSlljN  H£  T  A  (  I  I  )  *I>M(  1  1  ) 

R  f  A  £.  H «  L  *  L f  »  N X 
NR=0.u3b 

HttAWs?. 7* t | .-MM) /(NMF* i*i)R) 

DO  21  1=1 .*« 

fl(  1  A  <  I  1  =  2 , 7  »  f  1  , » »  (  I  )  )  /  f  'F  f  I )  *  ♦  I*  o  (  1 1  ) 

IMISO, 

TH2S0  , 

DO  4  J  =  1  »  ’* L 
Thjs7m1*Tm(J) 

NYl  s  J 

OM(J)slM< Jl/MO 

tf  £ 1M1  .U1  ,«0)  n»<  (  J  1  sCMn-TM?)  /MO 
If  (Imj  .i.t.moi  go  to  s 
Tm?s1h2f  T-( Ji 
SOM2=0, 

oo  ife  J=ifNru 

SOM|=fl. 

On  17  I  s i . n m 

sum)  s  s' ini  «af  TA(n/'’fT*»»t  (i.j) 
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mi*  nu  i*  ^S/**^*^ 

Won  39TI  tm*v..rr..  TOI»« 


1 


I 


16  Sum2ssU"?«0h(  J)  /( a  <  SO"! ) ) 

HI  UiHn 
f  Nf> 

SLtHHOuT  I*f  iNTf  KfN,T,L,nO,*,Mu.f),TT*RI«"t-,  IFL*G) 

COMM('N/Sf  t  I  /  N  R  l  ,  f  S 

CO*Mr>M/“AOS/F  ST  (<i,  111  ,DijT  (9,  1  t  )  ,  Rfd  7,1  1  )«TX(9,10)*RK(9«  101 
OI-lNSlUN  TS{ 10) .KSftO) 

RIAL  NKL.L.NU.KP.LAMHDa.N 
SS*(N/0,-iS)*«2 
KO*2.*3.,'AtS9/«L 
Nr  |  .N»R()«L 

«t  n^.MH  ,»m/(na* s*n> 

l  A”rl)A*t  . 

F  =  0. 

NC=1 70. 

I  c:  *  o 
2  f  Naf 

ic=  ic.i 

HsA«SljHT  (9,H0/M0)/(  1  ,»|  AMpOA) 
h  Dap*  r>/Nij 

fsn/(r(i*i  )*(SuHT(  i  ,♦(  i  .♦«c/»ri)*(  16.*»1  rA*A*L/(i,*5,iuiS9*Mn)))-i,j 

Lamhuasro*! *r/(2.*0 

lMlC.f.i.lM  M>  Hi  b 

lM(AF,S(»N-n/f ). i,i  no  to  ? 

S  T  Is  I  ,/(  I  ,*| 
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